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We report the discovery of a faint (Lx ~ 4 ± 1.5 x 10^''' erg s^^, 0.5-2 keV), 
out-flowing gaseous hot interstellar medium (ISM) in NGC 3379. This represents 
the lowest X-ray luminosity ever measured from a hot phase of the ISM in a 
nearby early type galaxy. The discovery of the hot ISM in a very deep Chandra 
observation was possible thanks to its unique spectral and spatial signatures, 
which distinguish it from the integrated stellar X-ray emission, responsible for 
most of the unresolved emission in the Chandra data. This hot component is 
found in a region of ~ 800 pc in radius at the center of the galaxy and has a total 
mass M~ 3±1 x 10^ M©. Independent theoretical prediction of the characteristics 
of an ISM in this galaxy, based on the intrinsic properties of NGC 3379, reproduce 
well the observed luminosity, temperature, and radial distribution and mass of 
the hot gas, and indicate that the gas is in an outflowing phase, predicted by 
models but not observed in any system so far. 

Subject headings: X-rays: ISM — galaxies: NGC 3379 — galaxies : elliptical and 
lenticular — X-rays : galaxies 



1. Introduction 

The nature of the X-ray emission of elliptical galaxies has been debated since their 
discovery as X-ray sources with the Einstein Observatory, which led to the early suggestion 
of ubiquitous hot X-ray emitting halos in hydrostatic equilibrium (Forman et al. 1985). 
These halos would originate from the stellar ejecta resulting from their normal evolution, 
accumulated during the lifetime of the galaxies, and retained by their deep dark matter 
potentials. However, from the analysis of these early observations and of ROSAT and ASCA 
data, a more complex picture of the X-ray emission emerged, suggesting that the non-nuclear 
X-ray emission of these galaxies consist of varying amounts of hot ISM, and of a baseline 
component related to the stellar population, either from the finale stages of stellar evolution 
(i.e. low-mass X-ray binaries, LMXBs, Trinchieri & Fabbiano 1985; Canizares et al. 1987; 
Kim, Fabbiano & Trinchieri 1992; Fabbiano et al. 2006; Matsushita et al. 1994) or possibly 
from coronal stellar emission (Pellegrini & Fabbiano 1994). Although controversial at the 
time, this picture has important implications for the evolutions of the stars and of the gaseous 
components of ellipticals, suggesting halo retention in the most massive systems, and partial 
or full winds in most galaxies (Ciotti et al, 1991; David et al. 1991; Pellegrini & Ciotti 1998). 

With Chandra for the first time the X-ray emission of nearby ellipticals can be imaged, 
and the existence of the different emission components can be pursued by direct observational 
means. These observations have led to the detection and study of LMXB populations (see 
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review in Fabbiano 2006), and the removal of their contamination from the diffuse emission 
component (e.g. NGC 1316, Kim & Fabbiano 2003). In some galaxies this diffuse emission 
is due to gaseous halos; however, in other galaxies, where the X-ray luminosity function 
can be probed to very low luminosities and the residual emission is itself of low luminosity, 
unresolved LMXBs and possibly stellar emission may account for a large amount or even the 
bulk of the X-ray luminosity (e.g.NGC821, Pellegrini et al 2007a,b; M32, Revnintsev et al. 
2007). 

We can now attempt to constrain the presence and properties of vestigial hot halos in 
these gas-poor galaxies. Since ellipticals are expected to host super-massive nuclear black 
holes (e.g., Richstone et al. 1998), which are however radiatively quiescent in most cases 
(Fabian & Canizares 1988, Pellegrini 2005, Ho 2008), these measurements provide useful 
constraints for models of halo evolution and nuclear feedback (e.g., Springel, Di Matteo 
& Hernquist 2005; see Pellegrini et al 2007b). The nearby, gas-poor, early-type galaxy 
NGC3379 provides an excellent target for these investigations. This is a very well studied 
prototypical elliptical, with Lb = 1.5 x 10^'^ Lb©. A series of very deep exposures were 
obtained for this galaxy by our team (PI: Fabbiano), as part of a very large Chandra program 
aimed at unraveling and studying the different X-ray emission components of ellipticals. 
Together with pre-existent archival data (David et al 2005), these observations resulted in 
an integrated exposure of 337 ks with ACIS S-3 (see Brassington et al. 2008 for details). In 
this paper we present our analysis of the diffuse/unresolved X-ray emission component of 
this galaxy. This analysis leads us to conclude that in addition to the unresolved emission 
of faint LMXBs and other stellar sources (Revnivtsev et al. 2008), a hot ISM is likely to be 
present in the central 800 pc of NGC3379. 

We discuss our data analysis in Section 2, and our interpretation of these results in 
Section 3. Our conclusions are summarized in Section 4. We assume a distance of 10.6 Mpc 
(Tonry et al. 2001), which gives a scale of 51 pc/". 



2. Analysis and Results 

We have used the same dataset that Brassington et al. (2008) used to determine the 
source list and compose the most complete X-ray catalog to-date in NGC 3379, down to 
(0.3-8.0) keV luminosities Lx ~ 10^^ erg s^^. All sources in the catalog have been excluded 
from the present analysis. We have used an aperture of 2" radius. On axis, this aperture 
encircles > 97/95% of the Point Spread Function (PSF), for energies 0.5/1.5 keV (see Figures 
4.6 and 4.20, and Table 4.2 in The Chandra Proposers' Observatory Guide, and our own 
simulations with ChaRT, see Appendix), and it degrades very slowly to encircle ~ 90% at 



-4- 



> 3' for the same energies (Fig. 4.13 in the The Chandra Proposers' Observatory Guide). 
As will become apparent in what follows, this work will concentrate on the central area 
NGC 3379, of ~ 80" radius; therefore, our data are not affected by the degradation of the 
Chandra point response function, and we do not need to consider a larger area for source 
exclusion. We discuss the possible contamination of the PSF in the Appendix. 

The size of the area considered also ensures a very homogeneous coverage from all five 
observations, since the region under study resides well within the area uniformly covered 
in all the Chandra exposures (roughly a 8' x 5' rectangle, see Fig. 1 in Brassington et al. 
2008). Therefore we do not need to resort to modelling of the exposure, and we can use a 
local estimate of the background from the same merged dataset. This procedure minimizes 
the uncertainties introduced by using a template or "blank sky observations" of different 
portions of the sky at different times, to estimate the field background. 

We have used two different merged datasets for our analysis. All five observations were 
used for the spatial analysis but only the four most recent observations were merged for the 
spectral analysis, given the significantly different response of the ACIS S3 CCD in the first 
observation obtained in 2001 (A02). Since this observation has the shortest exposure, not 
using it will not significantly degrade the statistics, while reducing the systematic error from 
the calibration uncertainties. The responses of the last four observations over the area of 
interest are virtually identical 



2.1. Radial Distributions of the X-ray Surface Brightness 

Fig. [T] shows azimuthally averaged radial profiles of the diffuse X-ray emission in the 
(0.3-2.0 keV) and (2.0-5.0) keV energy bands. The choice of the spectral ranges is motivated 
by our intention to discriminate between a harder band (2.0-5.0) keV, where all the emission, 
resolved and unresolved, is likely to arise from stellar sources, and a softer band (0.3-2.0) keV, 
to which a hot ISM may contribute. We have excluded all detected sources (as explained 
above), and centered the concentric annuli on the nucleus, identified as the optical/IR center 
and coinciding with source # 81 in Brassington et al. (2008). Note that as a consequence 
of this, and also due to source crowding at the center, the X-ray profiles cannot probe the 
very central region within r=2-3". The raw profiles shown in Fig. [1] flatten at a radius of 
> 80", indicating that the field background dominates. We can therefore use the 100" — 140" 
region to estimate the local background and subtract it from the emission, to produce the 
net profiles shown in Fig. [T]-right in the same energy bands. Note that the soft profile is 
more extended and more centrally peaked relative to the hard profile. 
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In Fig [2] we divide the full energy range considered in three, and show the azimuthally 
averaged net profiles in the (0.3-0.7) keV, (0.7-1.5) keV and (1.5-5.0) keV bands. This choice 
of energy boundaries was motivated by the results of the spectral analysis of the diffuse 
emission (see next section, FigE]), to maximize the contribution of the two optically thin 
thermal plasma component identified in the X-ray spectrum, with temperatures of ~ 0.3 
keV and at ~ 1 keV, respectively. It is clear from Fig. [2] that while both the very soft and 
the harder component follow a similar radial distribution, the soft (0.7-1.5 keV) component 
has a distinct excess for galactocentric radii r< 20". 

In the same Fig. [2] we plot the radial profiles in the optical and in the near IR K-band, 
from the HST F814 filter image and the 2MASS k-band image, obtained from their respective 
archives (see also Cappellari et al. 2006). These profiles are also azimuthally averaged, but 
do not cover the full 360° azimuthal plane of the galaxy, except for the very center (10" 
and 30" radius for HST and 2MASS respectively) since the two images only partially cover 
the galaxy. We also add the Sauron I-band data profile from Cappellari et al. (2006) that 
cover instead the full extent of the galaxy. The profiles are normalized to the X-ray data by 
rescaling them for radii r> 100". 

With the exception of the very central region (< 3" radius), the optical-IR profiles follow 
one another, and trace well the distribution of the very soft and hard X-ray profiles. On 
the contrary, the (0.7-1.5) keV profile is more centrally peaked and only becomes consistent 
with the shape of the optical-IR ones at r> 15". 



2.2. Spectral Analysis 

We have analyzed the spectral data extracted in the 2" — 15", 2"— 30" and 2"— 45" regions 
from the merged dataset from the four more recent observations. Most of the results discussed 
here come from the 2" — 30" region, which is a reasonable compromise that includes a large 
fraction of the diffuse emission at high significance, but all regions studied give a consistent 
picture. To increase the statistical significance of each spectral bin, we have binned the data 
to obtain a minimum significance of la or better in the net data. The spectral fits are done 
assuming the line-of-sight galactic N/^ of 2.7x10^*^ cm~^, the abundance tables of Wilms et 
al. (2000), the APEC model to account from the optically thin thermal emission of a plasma, 
with abundances fixed at the 100% value (different values do not change the final results), 
and a Bremsstrahlung to account for the X-ray emission of unresolved fainter LMXBs and 
stellar X-ray sources. 

We find that a single APEC component, with kT~0.3 keV, and a 7 keV Bremsstrahlung 
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model are able to approximate the data and give an acceptable value for the reduced x^ ~ 
1. However, inspection of the residuals shows a noticeable excess at around 1 keV (with a 
peak at about 5 a) and a less significant deficit at 0.6-0.8 keV, as shown in Fig. [31 The 
addition of a second APEC component reduces the x^ value and most of all eliminates 
the excess, suggesting two separate components at 0.3 keV and 1.0 keV (Fig. [3], bottom). 
Fig. m shows the 68%, 90% and 99% confidence regions for the two interesting parameters, 
the temperatures of the very soft and soft components. Fig. [5] shows the best-fit unfolded 
spectral components, from which it becomes evident that the ~ 1 keV component is clearly 
dominant over the ~ 0.3 keV component in the 0.7-1.5 keV band. 

We tested whether the existing data support a requirement of a second soft plasma 
component in the spectrum, by computing the Bayesian Information Criterion (BIG), dis- 
cussed by Schwarz (1978; an astronomical introduction to it can be found in Liddle 2004, 
and an application in Tajer et al. 2007) for both models. This quantity approximates the 
Bayes factor (Jeffreys 1961; Kass & Raftery 1995), which gives the posterior odds of one 
model against another, presuming that the models are equally favored prior to the data 
fitting (Liddle 2004). The difference ABIC can be used to evaluate the relative value of the 
two models. We find ABIC of 5-6 (depending on the region used), which can be used as a 
strong evidence (see Jeffreys 1961; Liddle 2004) that we are justified in introducing the 1 
keV component in the spectral modelling. 

We have investigated further the spatial distribution of the different spectral components 
by comparing the spectra from the 2" — 15" and 20" — 30" regions, both fitted with a single 
APEC spectrum at 0.3 keV (Fig. [6]). It is clear that there is no residual emission at 1 keV 
in the spectrum from the outer region; this suggests that the additional "soft" component 
is more concentrated towards the center, consistent with the radial profile comparison of 
Fig.El 

From the spectral results we can estimate the luminosity of each individual component, 
summarized in Table 1. The measured unabsorbed luminosity of the total unresolved emis- 
sion, within the 30" region, is Lx ~ 4 ± 0.4 x 10^^ erg s~^ in the 0.3-10 keV range, equally 
divided below and above 2 kev. The hard band luminosity is totally accounted for by the 
Bremsstrahlung component, which alone gives a luminosity L^. ~ 1 ± 0.1 x 10'^^ erg s~^ in 
the 0.5-2.0 keV band. The 1 keV component contributes for Lx ~ 41;q'2 x 10^^ erg s~^ in this 
band, and the residual Lx ~ 6lg'4 x 10^'^ erg s~^ is due to the 0.3 keV emission. Errors are 
estimated from the best fit values (68% confidence values, from XSPEC). We have estimated 
a correction factor to recover the area lost due to masking the point sources (see Table 1), 
and a second one to account for the fraction of emission not included in the area used in 
the spectral analysis (given in Table 2). These can be used to estimate the total luminosity 
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of the very soft and hard components, which we hst in Table 2. The component at 1 keV 
appears to contribute only in the central region, for r < 30". 



3. Summary of results and comparison with the literature 

We have analyzed the spectral and spatial distribution of the diffuse X-ray emission 
in NGC 3379 with deep Chandra data, cleaned of all detected sources down to a threshold 
luminosity L^^ ~ 10^^ erg s~^ (see Brassington et al. 2008). The residual emission can 
be detected only out to r ~ 1', corresponding to 3 kpc, well within the D25, measured at 
5.4' X 4.8' (from NED, see also Cappellari et al. 2006). 

The spectral distribution of this unresolved emission in the 30" region used for the 
spectral analysis suggests the presence of 3 distinct components. The dominant contribution 
in the broad 0.3-10 keV band can be modeled with a Bremsstrahlung emission at ~ 7 keV 
and is likely to arise from lower luminosity LMXBs and other stellar components that can 
not be detected individually in the present data, as also suggested by its spatial distribution, 
which follows that of the optical and IR stellar light (Fig. [2]). This component accounts for 
all of the emission above 2 keV and about 1/2 of the emission below 2 keV, and has an 
estimated total L^ ~ 5 x 10^^ erg s"^ (0.5-10 keV, see Table 2). 

We also find two additional optically thin thermal components, with best fit tempera- 
tures of ~ 0.3 keV and ~ 1 keV, contributing roughly 2/3 and 1/3 of the 0.5-2.0 keV residual 
luminosity, respectively (see Table 1), in the 30" radius region. This result is supported by 
the combined evidence of excess emission at ~ 1 keV over the dominant 0.3 kev plasma 
in the spectral data (Fig. [3]), and of a different spatial distribution of the two components, 
when the data below and above 0.7 keV are selected (Fig. [2]). 

On average, the 0.3 keV component has a spatial distribution consistent with that of 
the lower luminosity LMXBs and coronally active stars (traced by the emission above 1.5 
keV) and of the stars (traced by the I-band and K-band emission), while the 1 keV emission 
is more centrally peaked, and is prominent out to r ~ 15" — 20". 

Two previous detailed studies of the diffuse emission of NGC 3379 based on Chandra 
data can be found in the literature. Based on the first 2001 short observation David et al. 
(2005) conclude that the diffuse emission is dominated by undetected faint point sources, 
but with a 10% contribution from a gaseous component, with kT~ 0.6 keV, in the central 
770 pc. This hot ISM would have a luminosity of ~ 9 x 10^^ erg s~^ (0.3-10 keV) and a 
mass of ~ 5 x 10^ Mq. Our estimate of the luminosity of the optically thin thermal diffuse 
emission in the same region is consistent with that of David et al (2005), but we are able 



to distinguish two separate components of this emission. NGC 3379 is also included in the 
sample examined by Fukazawa et al. (2006), who also report a detection of a component at 
0.5 keV in a small region of the galaxy. 

More recently, Revnivtsev et al. (2008) have suggested a different interpretation for 
the unresolved emission: using the full dataset as in the present work, they derive a radial 
profile of the 0.5-2.0 keV emission, cleaned of detected sources, which they compare with the 
K-band profile. Based on this comparison they conclude that all of the unresolved emission 
observed in NGC 3379 can be attributed to the emission of stellar sources, with no evidence 
of a truly diffuse ISM. Our net radial profile of Fig. [1] is entirely consistent with that of 
Revnivtsev et al. (2008), giving us confidence that slight differences in the data analysis, in 
the subtraction of detected sources obtained with different algorithms or in the choice of the 
background level, are not significant sources of discrepancies. The main difference between 
this work and Revnivtsev et al. (2008) therefore lies in our spectral and spatial analysis, 
that enables us to distinguish two separate spectral components of the soft diffuse emission, 
with different spatial properties, which contribute equal amounts to the total in the inner 
15" radius. 



4. Discussion 

Our analysis suggests the presence of three components of the diffuse emission of NGC 
3379. As summarized above, the hard and very soft emission are spatially consistent with 
a stellar origin; a detailed discussion of this integrated X-ray emission from stellar sources 
(LMXBs, cataclysmic variables, active stars and stellar coronae) is given by Revnivtsev et 
al (2008). This stellar emission has been invoked to explain the diffuse X-ray emission of the 
Galactic Ridge (Revnivtsev et al 2006, Sazonov et al. 2006), the bulge of M31 (Bogdan & 
Gilfanov 2008), M32 (Revnivtsev et al. 2007) and NGC821 (Pellegrini et al 2007a). 

In what follows we will concentrate on the emission of the (0.7-1.5) keV band. This 
component has a different spatial distribution, more concentrated in the center and more 
centrally peaked than the stellar emission (see Fig. [2]), indicating a different origin. We 
suggest that this component is likely to originate from a ~1 keV hot plasma, and discuss it 
in the context of an appropriate model for a galactic flow. 
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4.1. Evidence for hot ISM in NGC 3379 

The interpretation of tlie 0.7-1.5 keV emission detected in NGC 3379 as due to a different 
component than the softer and harder ones opens the possibility that this is the long sought 
for evidence of a hot ISM in this galaxy (e.g., Ciotti et al. 1991, David et al. 2005). However, 
as shown by Fig. \5[ all three components contribute to the 0.7-1.5 keV emission, while both 
the softer, and in particular the harder, energy ranges are less contaminated by the other 
components. The spectral results can be reliably used to estimate the contribution of each 
component to the total luminosity. To estimate the spatial distribution of each separately, 
we have to resort to a few reasonable assumptions, which we have already discussed and 
that we summarize below: a) the profiles in both the very soft and hard energy bands are 
very similar in shape, therefore we can treat them as having a common origin with different 
rescaling factors; b) they are also consistent with the I-band and K- band profiles, again with 
possibly different rescaling along the full spatial extent probed; and c) the 0.7-1.5 keV profile 
is also consistent in shape with all above profiles for r> 15", where there is no longer any 
evidence of the emission due to the 1 keV component (Fig. [H]). We can therefore reasonably 
conclude that the 0.7-1.5 keV emission at r> 15" traces the soft and hard components, and 
that an extrapolation at smaller radii could be used to give a reasonable estimate of the local 
contribution from the stellar component at different radii. We have used the K-band data 
normalized to the 0.7-1.5 keV emission at r> 15" as the best proxy to extrapolate the profile 
at smaller radii, since we can use it as a model with virtually no errors. While this relies on 
the assumption that the stellar light is a good tracer of the X-ray binary population, with no 
local, though small, variations, the observed agreement with the X-ray profiles of the residual 
emission at softer and harder bands (Fig. [2]) supports this assumption. We therefore derive 
the radial distribution of the ~ 1 keV component as the excess over our model (Fig. [7]) 
which we can now compare to the expected emission from a hot ISM, derived from the 
hydrodynamical simulations described below. Note that the X-ray luminosity derived from 
integrating the radial distribution of Fig. [7]isLa;~2.5±lx 10^"^ erg s~^, consistent with 
that derived from the spectral analysis. 

From the spectral results and the assumption of spherical symmetry, we can estimate 
the average density and total mass of gas responsible for the emission measured. We derive 
a very low density of n^ ~ 0.007 cm~^ and a total mass of gas M~ 3 x 10^ Mq, in a region 
of ~ 770 pc in radius at the center. 
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4.2. The hot gas from hydrodynamical simulations 

The round optical shape (morphological type El), its optical colors, typical of an old 
population (Terlevich & Forbes 2002), and the lack of signs of current or past interactions 
(Schweizer & Seitzer 1992) allow us to assume that a simple "passively evolving" stellar pop- 
ulation feeds a hot gas flow, whose evolution can be reasonably studied with spherically sym- 
metric hydrodynamical simulations. The numerical code used to solve the time-dependent 
equations of hydrodynamics with source terms is described in details in Ciotti et al. (1991). 
We adopt a central grid spacing of 20 pc, which provides a good match to the scale of the 
observations in the inner regions; the simulations do not cover feedback effects and the flow 
evolution close to the central massive black hole (Shapiro et al. 2006). 

Given the high quality data available for this galaxy, it is possible to build a mass model 
tailored specifically onto NGC 3379, and consider the best recipes available for the mass and 
energy source terms. These are briefly described below. 



4.3. The mass model 

The optical profile is well described by a de Vaucouleurs (1948) law over a span of 10 
magnitudes (Capaccioli et al. 1990; Peletier et al. 1990), therefore the simulations use the 
stellar mass density profile given by the Hernquist (1990) distribution, that is a very good 
approximation of the de Vaucouleurs law and has the advantage that its dynamical properties 
can be expressed analytically. We also impose the observed B-band luminosity of 1.5x10^° 
Lb0, the central stellar velocity dispersion of 230 km s~^ (from the detailed modeling of 
ground-based and HST spectroscopy, Shapiro et al. 2006) and a total stellar mass M* = 
lO^^M©, obtained from dynamical and stellar population synthesis studies (Saglia, Bertin & 
Stiavelli 1992, Gerhard et al. 2001, Napolitano et al. 2005, Cappellari et al. 2006, Douglas 
et al. 2007). The resulting shape of the stellar profile, when projected, is close to the optical 
profiles in Fig. [2l 

The radial density distribution of the dark haloes of ellipticals is not well constrained by 
observations; theoretical arguments and high resolution numerical simulations of dissipation- 
less collapse produce a density distribution ex r~^ near the center (Dubinsky & Carlberg 
1991, Navarro, Frenk, & White 1996). The model galaxy is therefore a superposition of two 
Hernquist density distributions, one for the luminous matter and one for the dark matter. 

The total amount of mass in NGC 3379 has been calculated several times using the 
observed stellar velocity dispersion profile, extending out to ~ 1 — 2 effective radii (i?e, e.g., 
Saglia, Bertin & Stiavelli 1992, Kronawitter et al. 2000, Samurovic & Danziger 2005), and 
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dynamical tracers as planetary nebulae or globular clusters extending out to larger radii 
(Romanowsky et al. 2003, Dekel et al. 2005, Teodorescu et al. 2005; Pierce et al. 2006, 
Bergond et al. 2006, Douglas et al. 2007, De Lorenzi et al. 2008). Regardless of the details 
of the modelling, studies based on the observed stellar velocity dispersion profile invariably 
show very little dark matter within Re, or even no dark matter. Within 2-3 Re, the stellar 
component alone satisfies reasonably well the observations, and the amount of dark matter 
is constrained to be at most very modest (with a dark to stellar mass ratio of the order 
of Mh/M^ ~ 0.1, Kronawitter et al. 2000; Cappellari et al. 2006). The evidence for dark 
matter is at the largest scales (10-40 kpc), probed mostly by globular clusters (see also the 
analysis of an HI ring, Schneider 1985), with a global galactic value of M^/M^ = 1.1 — 1.5 
(Dekel et al. 2005), Mh/M^ ~ 2 (Puzia et al. 2004), up to Mh/M^ ~ 6 (Bergond et al. 
2006). 

In order to make the model mass distribution consistent with the observations, we also 
imposed the ratio M,/M. = 0.1 within one rR while on the global galactic =cale it ha= 

been varied instead between 1 and 6. 



4.4. Time evolving inputs 

The model does not assume a steady state configuration, therefore the input ingredients 
of the numerical simulations, which are the rates of stellar mass loss from the aging stellar 
population and the rate of SNIa heating, are both evolving with time. In the numerical 
code the exact mass return rate prescribed by the stellar evolution theory is used (Ciotti et 
al. 1991), updated to take into account the latest stellar population synthesis models and 
more recent estimates for the mass ejection from stars as a function of their mass (Maraston 
2005). This gives a present epoch (~ 10 Gyr) mass return rate of ~ O.3M0 yr~^ for the 
whole galaxy. 

The SNIa heating rate is parameterized as LsNit) = Esn RsN^t) Lb, where Esn is the 
kinetic energy injected in the ISM by one SNIa, and the number of events as a function of 
time is RsN(t)LB oc t~*, where the slope s describes the unknown decay rate. RsN(t) is 
normalized to give the present SNIa's explosion rate in nearby E/SOs of Cappellaro et al. 
(1999), i.e., 0.16/ifo SNu (where 1 SNu = 1 SNIa per 100 yrs per 10^°Lb,q, h^ = HJ70), that 
has an associated uncertainty of ~ 30%. A detailed theoretical modeling of the evolution of 



^In the code the presence of dark matter contributes also to produce the observed central stellar velocity 
dispersion, as described in Pellegrini & Ciotti (2006), though this aspect is not very important here given 
the small amount of dark matter in the central regions. 
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the rate shows that, after the first 0.5-1 Gyr, it is well approximated by a power law with a 
slope s ~ 1.0 — 1.2 for double degenerates exploders, and a slope s ~ 1.5 — 1.7 for the single 
degenerate exploders (Greggio 2005). 

The thermalization of the stellar mass losses to the local "temperature" set by the stellar 
velocity dispersion is another heating source, though of lower relevance. Its contribution is 
determined at each radius using the velocity dispersion profile obtained by solving the Jeans 
equation for the two-component Hernquist model in the globally isotropic case (Pellegrini 
& Ciotti 2006). For realistic anisotropy distributions, the difference with the isotropic case 
is small; in addition, the heating due to thermalization of stellar motions is significantly 
smaller than SNIa heating. 



4.5. Results and comparison with observations 

With the SNIa's explosion rate from Cappellaro et al. (1999) and the adopted mass 
model, the resulting gas flow phase is of a global outflow, with velocities ranging from 10^ 
to a few 10^ cm s~^, going from the central region to the galaxy outskirts (e.g., subsonic 
out to r~l kpc and supersonic outwards). When varying the total amount of dark matter 
within the range M^/M^ = 1 — 6 allowed for by observations, this flow pattern keeps basically 
the same (provided that the SNIa's rate increases together with Mh/M^,, but still within its 
observational uncertainty). 

For the resulting gas flows, the X-ray luminosity in the energy bands of 0.3-2 keV and 2- 
8 keV, and the temperature weighted with the emission in the 0.3-8 keV band (< kT >), have 
been calculated. In order to do this, the values of the cooling function A(T, Z), as a function 
of temperature, have been obtained with the APEC code in XSPEC, for solar abundance, 
consistent with the estimate for the stellar population (Terlevich & Forbes 2002); a "best 
flt" A was then derived from these values and used in the Lx and < kT > calculation (that 
takes into account the whole density and temperature distributions over the computational 
grid) . 

After ~ 9 Gyr, the gas has a 0.3-2 keV luminosity of 2 x 10^^ erg s~^, in good agreement 
with the observed one. This is a reasonable time scale that corresponds to the estimated 
age of NGC 3379: Terlevich & Forbes (2002) give an age of 9.0l^;9 Gyr; more recently, a 
mean age in the range 8-15 Gyr has been conflrmed by Gregg et al. (2004) and Denicolo 
et al. (2005). The central temperature is ~ 10^ K out to a few hundred pc, and is steeply 
declining outwards; the central gas density is rie = 2 x 10~^ cm~^ and is also steeply declining 
outwards. The mass of the hot gas is ~ 3 — 4 x IO^Mq within a radius of 800 pc (in good 
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agreeinent with that obtained in Sect. 4.1). The emission weighted gas temperature within 
a radius of 15" is < kT >~ 0.8 keV, a value near the lower range of those allowed for by 
observations (see Fig. H]). The projection along the line of sight of the 0.3-2 keV emission 
is shown in Fig. [TJ for two models with a total Mh/M^, = 4, a SNIa's rate close to that of 
Cappellaro et al. (1999), and a slope for the SNIa's decay rate of s = 1.0 and s = 1.5. Both 
model profiles reproduce well the observed "gas profile" derived above. We note that a larger 
age for the gas flow would result in a lower stellar mass loss rate and lower gas content, with 
consequently lower ISM emission, below the observed one. 



5. Conclusions 

With the very deep Chandra observation of NGC 3379 we have scored two firsts: we have 
the first evidence of a very small amount of hot ISM in a nearby galaxy, down to the level of 
Lx ~ 4 X 10^^ erg s~^ (0.5-2.0 keV), and at the same time we have the first evidence that this 
gas is in an outflow phase, which has been predicted by models but as of yet has not been 
observed. As a third and equally important result, we find that detailed modelling based on 
the optical data (stellar mass, total mass and supernova rate) reproduces remarkably well 
the observed total luminosity, total mass and spatial distribution of the hot phase of the gas. 

After detection and removal of point sources down to a luminosity of Lx ~ 10^^ erg s~^ 
we have been able to study the unresolved emission in NGC 3379. The detailed spectral 
analysis has led us to the evidence of three separate components: a very soft (plasma at 
kT~0.3 keV) and a hard (Bremsstrahlung at kT=7 keV) components can be attributed to 
the integrated emission of unresolved sources, closely linked to the stellar population (binary 
systems and coronally active stars, as discussed by Revnivtsev et al. 2008). However we 
also detect a hot ISM in this galaxy, within the innermost ~ 800 pc, with a luminosity 
Lx ~ 4 ± 1.5 X 10^^ erg s~^ (0.5-2.0 keV). This is lowest detection of a hot gas in a nearby 
early-type galaxy. The total gas mass involved is also very small, of a few 10^ Mq. 

With the aid of hydrodynamical simulations and a galaxy model tailored specifically on 
NGC3379, we have reproduced both the luminosity and the average radial distribution of 
this hot gas component. For a SNIa's explosion rate consistent with current estimates, the 
gas is outflowing, even for the "maximum" dark matter model, where the amount of dark 
mass is the maximum allowed for by optical studies. 
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6. Appendix: Consideration upon the effects of the PSF 



At the suggestion of the referee, we have thoroughly checked that the resuhs obtained 
above, in particular the detection of the ~ 1 keV component, are not due to a faulty sub- 
traction of the detected sources. We have run several tests, which we illustrate below. 

1) We have first redone the spatial and spectral analysis assuming a circle of 2.5" and 3" 
radius for point source exclusion, with analogous results, albeith somewhat larger errors, due 
to the smaller statistics available. 

2) We have also simulated the Chandra^s PSF using ChaRT (Carter et al. 2003). Since we 
were interested in the wings of the PSF, we have run several simulations for 20,000 counts 
each and then merged the resulting images to increase the statistics at radii larger than 
1" We have used the average spectrum obtained in the inner 15", inclusive of all detected 
sources, as an input parameter in all simulations. We have then treated the results as we 
have the data of NGC 3379, namely we have produced radial profiles in the same energy 
bands and calculated the expected contribution from point sources in the regions of interest. 
We expect that 4% of the counts detected from a point sources within 2" from its peak 
falls in the 2" — 30" region, and 3.6% in the 2" — 15" region, in the 0.3-2.0 keV band (with 
values between 3% and 5% in the softer and harder of the bands considered). If we integrate 
the contribution due to the "wings" of all point sources detected in a radius=30", we can 
only account for 10-25% of the detected counts in the 0.3-0.7, and 0.7-1.5 keV bands in the 
unresolved component. This is a generous estimate of the contribution from the discrete 
population, since only sources at the center of the region will contribute the full amount, the 
wings of the others being only partially included in the region considered. 



3) We have used the current documents onlinci and directly contacted the experts on both 
the PSF and its wings (Dr. T Gaetz and D. Jerius) to better understand the effects of the 
wings of the point spread function. The most updated studies on the PSF wings are based 
on the ground XRCF and on-orbit Her X-1 data. Both data sets demonstrate the existence 
of the PSF wings, but the effects of this component decrease rapidly with decreasing energy. 
The ray trace model used to simulate the PSF (point 2 above) includes some of the effects 
of wings, and appears to be reliable at low energies (< 2 keV) and small radii, although 
it might be underestimating the wing effects at energies above 2.0-2.5 keV. However, we 
are mostly interested in the emission below ~ 1.5 keV and relatively small radii, where the 
model appears to be adequate. 



http : //cxc . harvard . edu/cal/Hrma/psf /index . html 
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4) We have re-analyzed the spectrum of the emission in the 2" — 15" regions excluding 
only the 5 brightest sources with a significantly smaller exclusion region of 1" radius. This 
exclusion is necessary to avoid that the point source population completely dominates the 
spectral results; on the other hand, if the sources were responsible for the 1 keV component, 
it should be significantly stronger when sources are included in the spectral analysis. The 
spectral results indicate that, while the net counts increase by a factor of 4, the 1 keV 
component only increases by a factor of 1.5 the original value, consistent with the ratio of 
the areas of the regions considered. 

All in all, we believe that the excess emission observed in the 0.7-1.5 keV band is a real 
effect, and is not an artifact related to the wings of the many individual sources detected 
in the area. This excess is relative to the emission in adjacent energy bands, and is seen at 
energies where the effects of PSF wings are significantly smaller than at higher energies. We 
also believe that there is excess emission over the scattering from the detected point sources, 
although we are fully aware of the fact that the exact evaluation of its strength has large 
uncertainties, both in the soft and in particular in the hard band. 

We thank Stefano Andreon for many interesting discussion on the statistical aspects of 
the data analysis, and Terry Gaetz and Diab Jerius for their help in understanding the PSF 
issues. We thank the referee for raising the issue of the effects of the PSF, that prompted 
us to better investigate the reliability of the results. The data analysis was supported 
by the CXC CIAO software and CALDB, and has made use of the SAOImage DS9 and 
funtools softwares, developed by the Smithsonian Astrophysical Observatory. We have used 
the NASA NED and ADS facilities, and have extracted archival data from the Chandra, 
HST and 2MASS archives. This work was supported by the Chandra GO grant G06-7079A 
(PI: Fabbiano) and sub-contract G06-7079B (PI: Kalogera). We acknowledge partial support 
from NASA contract NAS8-39073 (CXC); A. Zezas acknowledges support from NASA ITS A 
grant NAG5-13056; S. Pellegrini acknowledges partial financial support from the Italian 
Space Agency ASI (Agenzia Spaziale Italiana) through grant ASI-INAF 1/023/05/0. 
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NGC 3379 - raw data 
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Fig. 1. — Azimuthally averaged radial profiles of the unresolved emission in NGC 3379, in 
two broad energy bands. Left: raw data; right: background subtracted data 
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Fig. 2. — Azimuthally averaged net profile in different X-ray bands, compared to the optical 
and K-band profiles 
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Fig. 3. — Spectral data in the 30" region, fitted with a single temperature (0.3 keV) and 
a two temperature APEC model (top and bottom respectively, see text), plus a 7 keV 
bremsstrahlung. For displaying purposes, data are rebinned to 5a. 
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Fig. 4. — Confidence contour regions (at 68, 90 99% level for two interesting parameters) for 
the two temperatures of the APEC component (see text). 
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Fig. 5. — Region 30", unfolded spectrum: 0.3 + 1.0 keV APEC + 7 keV Bremsstrahlung, 
line-of-sight absorption. 
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data and folded model 
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Fig. 6. — Spectral data in the 0"-15" and 20"-30" region, fitted with a single temperature 
plasma model 



-25- 




IQi 
R(arcsec) 



10< 



Fig. 7. — The observed 0.3-2 keV emission from the "pure" 1 keV component, compared 
with the profiles resulting for two gas models, one with s = 1.0 (dashed line) and the other 
with s = 1.5 (solid line). See Sect. IHfor more details. 
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Table 1: Measured luminosities of the different components 



L:-b 



T b 



T b 



Region Net counts'^ Ijx" i^x" -i^x" 

0.3-10 keV 0.5-2.0 keV 2.0-10.0 keV 



B/A1/A2= 
0.5-2.0 keV 



Correction'^ 



2"-30" 

2"-15" 



1550±53 
978±35 



4.4x10^8 
2.9x10^^ 



2.0x10^^ 
1.3x10^^ 



2.0x10^8 
1.3 xlO^^ 



1/0.6/0.35 1.15 

0.8/0.27/0.28 1.34 



"Net counts used in the spectral fit 
''Total luminosity in the given band in erg s~^ 

=Lx in the 0.5-2 keV band for the Brcmsstrahlung / APEC @ 0.3 keV/ APEC @ 1 keV components, in units 
of 10^^ erg s^i 
''Correction factor to take into account the area lost due to the excised sources 



Table 2: Total luminosities of the different components, after correcting for area lost and 
total extent 



Component 



Extension Correction^ 
farcsec) 



0.5-2.0 



Lx 
2-10 



Very soft (0.3 keV) 


95 


Soft (1 keV) 


20 


Hard (7 keV) 


50 



1.60 


1.1x10^^ 


2.3x10^5 


1.00 


4.0x10^^ 


5.0x10=^6 


1.35 


1.6x10^^ 


3.0x10^8 



"Correction factor to extrapolate from the spectral extraction region (30", see Tab. 1) to the total source 
counts 



